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ABSTRACT 

The protostellar collapse of a molecular cloud core is usually accompanied by outflow 
phenomena. The latter are thought to be driven by magnetorotational processes from 
the central parts of the protostellar disc. While several 3D AMR/nested grid stud- 
ies of outflow phenomena in collapsing magnetically supercritical dense cores have 
been reported in the literature, so far no such simulation has been performed using 
the Smoothed Particle Hydrodynamics (SPH) method. This is mainly due to intrin- 
sic numerical difficulties in handling magnetohydrodynamics within SPH, which only 
recently were partly resolved. In this work, we use an approach where we evolve the 
magnetic field via the induction equation, augmented with stability correction and 
divergence cleaning schemes. We consider the collapse of a rotating core of one solar 
mass, threaded by a weak magnetic field initially parallel to the rotation axis so that 
the core is magnetically supercritical. We show, that Smoothed Particle Magnetohy- 
drodynamics (SPMHD) is able to handle the magnetorotational processes connected 
with outflow phenomena, and to produce meaningful results which are in good agree- 
ment with findings reported in the literature. Especially, our numerical scheme allows 
for a quantitative analysis of the evolution of the ratio of the toroidal to the poloidal 
magnetic field, which we performed in this work. 
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1 INTRODUCTION 

Observations in star forming regions have revealed that the 
protostellar collapse of a molecular cloud core is usually ac- 



companied by outflow phenomena (e. g. Wu et al. 2004 



Bally|2007|), which were already discovered back in the 1970s 



(jZuckerman Palmer|1975| [Kwan Scoville|1976| |Zucker- 



man et al.|1976 ) . Usually, protostellar outflows are classified 
into two types, optical jets or molecular outflows. The lat- 
ter typically exhibit slow velocities and wide opening angles, 
and are observationally identified by line emission from their 
CO molecules. The former, having narrow opening angles 
and high velocities, are observed optically. Both phenomena 
may occur together, where usually the high-velocity jet is 
enclosed in a low-velocity molecular outflow. The physical 
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mechanisms that drive outflows are still not fully under- 
stood. The adiabatic (or 'first') core, the circumstellar disc 
around the protostar and the protostar itself are possible 
origins of either jet or molecular outflow or both. First an- 



alytical work by Blandford & Payne ( 1982 ) and Pudritz & 
Norman (1983) suggested rapid rotating, magnetised pro- 



tostellar discs as possible origins of these outflows, since 
other sources like pressure due to either the gas or radiation, 
have been ruled out since those mechanisms fail to provide 



sufficient energy and momentum. Later, Lynden-Bell ( 1996 



2003) introduced the idea of a 'magnetic tower', an outflow 
driven by magnetic pressure originating in toroidal magnetic 
field components. 

From a numerical point of view, the simulation by 
Tomisaka ( 1998 ) was one of the first to produce a magneti- 



cally driven outflow during the collapse of a magnetised and 
rotating molecular cloud core. There it was shown that the 
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outflow is launched by dynamically building up of a toroidal 
magnetic field. Other simulations concentrating on outflow 
phenomena were performed by |Tom isaka (2000 2002), and 
|Banerjee &; Pudritz| ( |2006| ) who followed the evolution up to 
protostar formation using an ideal MHD approximation and 
appropriate cooling models. These works in general agree in 
finding two distinct flows emerging from different regions, 
namely the low- velocity flow from the first, the high- velocity 
flow from the second core which, respectively, can be iden- 
tified as molecular outflow and optical jet. Another work by 



Hennebelle &; Fromang (2008) used also ideal MHD and a 
barotropic equation of state and found a magnetic outflow 
best described by a magnetic tower for a model with low 
magnetic fields strength, and a morphological different flow 
for a high field strength model with characteristics of the 
disc-wind model proposed by Blandford & Payne (1982). 



However, their barotropic equation of state only had a pre- 
scription for the formation of the first core, and so they 
were unable to model the evolution beyond this point in 
the collapse. A similar work by Machida et al. (2005}, how- 
ever starting from a filamentary cloud with initial density 
perturbations but also using a simple barotropic equation 
of state, also reported a low- velocity outflow from the first 



Machida et al. (2006) used a more sophisticated equa- 



tion of state that approximates the thermal evolution better, 
in particular it allows to follow the collapse until second core 
formation. Furthermore, they used ideal and resistive MHD 
models, and found a high-velocity jet ejected from the sec- 
ond core in both models. In a following and more detailed 



study, Machida et al. (2008) also stressed the importance 



of non-ideal MHD in high density regions. They considered 
models with different initial angular velocities and also found 
a low- velocity flow from the first core and high- velocity flow 
from the protostar. However, these authors suggest the for- 
mer flow being driven by disc-wind mechanisms and the 
latter by magnetic pressure, conversely to other authors. 



Machida et al. ( 2010| used again a simple equation of state 



Machida et al. (2005), but in combination with resis- 



tive MHD, and applied this scheme to a spherical cloud with 
initial density perturbations. They concentrated, however, 
on the formation of the circumstellar disc. They found an 
outflow driven by the circumstellar disc which, however, is 
weakened in the later main accretion phase since the amount 
of inf ailing gas decreases over time. Finally, the outflow is 
found to disappear. In a more recent study, which considered 
the combined effects of MHD and radiation and follows the 



collapse until first core formation, Tomida et al. ( 2010 ) found 



a larger size of the first core and higher entropy in its outer 
layers, compared to MHD simulations using a barotropic ap- 
proximation. Furthermore, they find a two-component out- 
flow driven by different mechanisms, namely an inner mag- 
netic pressure mode and an outer magneto-centrifugal mode. 
In this work, we perform SPMHD simulations using a simple 
set-up with initial parameters comparable to those chosen 



by Hennebelle Sz Fromang| ( |2008| and we concentrate on a 
weak field scenario. We show, that our results agree qualita- 
tively well to those obtained by these authors. In particular, 



Hennebelle & Fromang ( 2008 ) also suggest toroidal magnetic 



pressure as driver of the low- velocity outflow, in agreement 
with our results. 




Figure 1. Magnetic field strength B as a function of the cen- 
tral number density n c in the core. Each data point represents 
a spatial average from simulation snapshot which were taken in 
intervals of 0.01 %. The solid line was obtained from the power 
law B oc n K with n = 0.58. 



2 NUMERICAL METHOD AND INITIAL 
CONDITIONS 

We perform three dimensional simulations using the 
Smoothed Particle Hydrodynamics (SPH) code GADGET 
( |Springel et al.|2001| |Springel 



2005), which was recently ex- 



tended by |Dolag &; Stasyszyn ( 2009 ) to allow for the treat- 
ment of ideal magnet ohydrodynamics. We refer to the latter 
publication and references therein for a detailed discussion of 
the algorithms and numerical tests. This implementation of 
Smoothed Particle Magnetohydrodynamics (SPMHD) has 
been successfully applied in extra-galactic (iBonafede et al. 



2011) and galactic astrophysics ( |Kotarba et a l. 2009, 2010 



2011), as well as to collapse and fragmentation problems in 
star formation ( Biirzle et al.||20lT ). All parameters used to 
control the SPMHD scheme have the same values as adopted 
(2011), except that a = 1 in this work. 
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In this work, we use initial conditions reproduced from 



those given in the work by Hennebelle & Fromang (2008) 



The initial cloud mass was chosen to be Mo = IMq, the 
radius was set to Ro = 4.6 • 10 16 cm w 0.015 pc. This choice 
corresponds to an initial density of po = 4.87 • 10 -18 gem -3 
and a free-fall time ts ~ 30000 yr. The spherical cloud is em- 
bedded into a low density medium with 5 = po/pext = 100. 
This medium, however, is not in pressure equilibrium with 
the cloud (Po/^Pext = 100) thus allowing the outer cloud 
layers to expand into the medium, as in Hennebel le &; Fro^] 
mang ( 2008| . The whole system is placed into a cube with 
side length of 4Ro and periodic boundary conditions. 

The thermodynamics is controlled by a barotropic equa- 
tion of state given by 



2 

c s ,oP 



1 + 



VPcrit / 



4/3 



1/2 



(1) 



where c s ,o = 0.2kms _1 is the speed of sound, result- 
ing from a temperature T — 11 K and a mean molecular 
weight of U- m oi = 2.1. The critical density is chosen to be 
y9 cr it = 10 -13 gem -3 , leading, as a consequence of the longer 
remaining of the gas in the quasi-isothermal regime, to a 
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very rapid growth of the central density. Thus it is difficult 
to follow the evolution for a long period of time, especially 
the second collapse is not within the range accessible in this 
work. Usually this problem can be circumvented by using 



sink particles ( 


Bate et al.||1995 Jappsen et al.||2005| Feder- 


rath et al. 2010 


) , but since we investigate outflow phenomena 



launched in the proximity of the 'protostar' we avoid using 
them here to reduce perturbations to the flow behaviour as 
far as possible. 

We define the ratio of the mass-to-flux ratio to the crit- 
ical mass-to-flux ratio as \i — (M/$)/(M/$) cr i t , where <3> is 
the magnetic flux. For a spherical configura tion, the critical 
value is given by (M/$) C rit ~ 0.125 G~ 1/2 (|Mouschovias h 



Spitzer 1976). Choosing \i = 20 the initial magnetic field 
strength is given by 30.7 |^G, and the plasma beta is given 
by /3piasma = 56. The Alfven speed is va = B / ^Anpv = 
3.9 x 10" 2 kms _1 . 

The initial angular velocity is = 4.3 • 10 _13 s _1 for 
the cloud which is in solid body rotation, corresponding to 
£ff = 0.4. Using these choices, we get for the energy ratios 

Q^therm — ^therm/ | -E'grav | — 0.37 and /3 ro t = E TO t / \ E gr &v \ = 

0.045. Machida et al. (2008) characterize their models by a 

Qo / y/AirGpo resulting for our 
value of (jj = 0.2, thus com- 



dimensionless parameter us 
choice of initial conditions in 
parable with the rapid rotating model considered by these 
authors. 

The Jeans condition for SPH given by Bate & Burkert 



(1997) requires that the mass of the cloud needs to be suffi- 
ciently resolved, and that the particle softening length needs 
to be equal to the smoothing length, in order to avoid numer- 
ical artefacts during cloud collapse. They proved, that the 
minimum number of particles required to resolve the cloud 
mass is given by iV m in = 2MoAf ne igh/Mj, where Mj is the 
Jeans mass. There exist several formulations of the Jeans 
mass which differ, as was pointed out by Nelson (2006), by 
numerical prefactors. We use the most conservative formu- 
lation that gives the lowest value for the Jeans mass: 



Mj 



(2) 



1/2 /5\ 3/2 cj 
W G 3 /V /2 " 

For cl/p~ 1/2 at p crit we get, using eq. |l]), 2 3/4 c^op~^ 2 and 
thus a value of N m i n ~ 48000 for N ne i g h = 64. Employing 
N = 4 • 10 6 particles within the cloud (6 • 10 5 in the am- 
bient medium), we have N ~ 83iV m in and thus meet the 
requirements imposed by the Jeans condition. Additionally, 
we would like to emphasize that we also did simulations with 
~ 1/5 of the particles and found qualitatively the same be- 
haviour as described here. 



3 RESULTS & DISCUSSION 
3.1 Scaling of the magnetic field 

Regarding the evolution of the magnetic field, we show in 
Fig. |l]) the dependency of the magnetic field B on the 
core density n c . Assuming a scaling relation parametrized 
as B oc n K , as is usually done in the literature (e. g. iHeilesI 



&; Crutcher 2005), we performed a fit which gave a result 
of k = 0.58. For this fit we considered only the region of 
the highest density n m ax within the collapsing cloud. In- 
creasing the considered density range used for averaging up 



to rimax/2 ^ n ^ rimax, as in Banerjee & Pudritz (2006), 
we see the same trend only slightly steeper with k = 0.64 



(not shown). This is in good agreement with Banerjee & 
|Pudritz| ( [2006t and |Price fe Batej (2007) , who find a v alue of 
■ 0.6, but in disagreement with other studies (e. g. Desch 



&; Mouschovias|2001|p et al.|2004||Burzle et al.|2011| ). How- 
ever, the scaling relation depends strongly on the geometry 
of the magnetic field. 

3.2 First core formation and outflow 

We define the adiabatic, or 'first', core as the space occupied 
by the material exceeding the critical density p C rit, that is 
where the gas becomes optically thick. We find that the 
first core has formed at a time of 1.1 1& or 33000 years. At 
this time the central density has reached a value of n max — 
10 11 cm -3 , and the mass of the first core is M « 0.02 M . 
Due to the fast rotation the first core has an oblate structure, 
with an equatorial radius of r eq = 43.5 AU in the x — y plane, 
and a polar radius of r po i = 9.5 AU parallel to the z-axis at 
formation time. Defining the oblateness as usual, we find a 
value of £ = 1 — r P oi/r e q = 0.78. 

Only shortly after formation of the first core, at ~ 
1.11 tff, a slow bipolar outflow with peak velocities of ~ 
1.9 km s -1 is ejected from the disc. This can be seen in 
Fig. |2]), which shows the normalized density near the x- 
z plane for different timesteps, and the vectors show the 
velocity field (arrows have equal magnitude for better visi- 
bility). This outflow is accompanied with the build-up of a 
strong toroidal magnetic field, as shown in Fig. (|3|, where 
the ratio B t or/B po i > 1 in the outflow regions. This kind 
of low- velocity outflow, which is frequently found in other 
simulations (e. g. [T omisaka 1998; Ba nerjee &; Pudritz| 2006 



Hennebelle fe Fromang||2008| |Machida et al.||2008|), is usu 
ally considered as magnetic tower (Lynden-Bell 1996, 2003 



Kato et al. 2004). In a magnetic tower, the toroidal com- 
ponent of the magnetic field is continuously generated by 
the rotating disc and pushes material into the surround- 
ing medium by means of magnetic pressure. Additionally, 
this can be proven by considering the ratio of thermal to 
magnetic pressure which, in the outflow region, is given by 
/^plasma ~ 0.1 indicating magnetic pressure dominating ther- 
mal pressure there. This magnetic tower flow creates a torus- 
like structure that further expands into the surrounding gas, 
also visible in Fig. |2]). 

We note that our results suggesting toroidal magnetic 
pressure as driving mechanism are in good agreement with 
results presented by |Banerjee & Pudritz| ( |2006| and |Hen 



nebelle & Fromang (2008), but, however, in disagreement 
to the findings by|Machida et al. ( 2008 ) . The latter authors 



attribute the emergence of the slow outflow to magneto- 
rotational effects rather than magnetic pressure. On the 
other hand, they had included the effects of non-ideal MHD 
(ambipolar diffusion and resistivity) which are absent in this 
work and in|Banerjee &; Pudritz (2006) and Hennebelle & 



Fromang (2008) which might have an important influence 
on the magnetic field structure. 

Later in our simulation, at 1.17 tff, also a fast outflow, 
with peak velocities of ~ 28kms _1 , emerges from the cen- 
tral parts of the first core, visible in the last panel of Fig. 
|2]) an d Fig. (|3|), respectively. A similar structure is also vis- 
ible in Hennebelle & Fromang (2008) in the weak field case, 
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v 



Figure 2. The figure shows a sequence of panels with increasing time, starting at 1.10 % where the first core forms. The subsequent 
snapshots were taken at 1.12 %, 1.14 %, 1.15 %, 1.16 % and 1.17 %. Shown is the normalized projection of the density, i. e. an integration 
along the line of sight divided by the average value, performed in a thin slice around the x-z plane. The arrows which, for better visibility, 
were set to equal magnitude, indicate the velocity. 




PQ 

V 



Figure 3. The same sequence of panels as in Fig. |2j is shown. Here, the normalized projection of the ratio of the toroidal to the poloidal 
magnetic field component Btor/^pob is shown. The vectors -Btor and B po \ were inferred from B = (B x , B y , B z ), which is known at 
every particle position, by expressing the projection of B onto cylindrical unit vectors in rectangular coordinates. The integration was 
again performed in a thin slice around the x-z plane. 



probably also at the same simulation time. However, since 
we use a comparable set-up as these authors, the same re- 
strictions regarding the thermal structure of the protostar 
also apply to our work. So, although the velocities found 
here are comparable to those obtained in other work (e. g. 
Banerjee fe Pudritz||2006| |Machida et al.|2008) , an identifi- 
cation of this fast flow component with optical jets seen in 
observations, which believed to be launched during second 



collapse, seems to be doubtful at this stage. Certainly, fur- 
ther investigations considering a more correct treatment of 
the thermodynamics within molecular cloud core are neces- 
sary. 

Finally, in Fig. ^ we show a sequence of spatial aver- 
ages of the ratio Bt OY /B po i from 1.10 t& to 1.16 tff, calculated 
in the central parts of the cloud core. Here we see the gen- 
eration of the toroidal part of the magnetic field, as well as 
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Figure 4. Shown is the spacially averaged ratio Bt OT /B po \, plot- 
ted against the distance r from the cloud centre for different sim- 
ulation times. As in Fig. J3I, Btor and B po \ were obtained by 
expressing the projection of B onto cylindrical unit vectors in 
cartesian coordinates. The curves show generation of the mag- 
netic field components and their outward propagation. 



its outward propagation with time. However, at 1.15 it 
is also visible that the toroidal part is re-generated in the 
very central parts of the disc. Thus, this proves the connec- 
tion of the outflow with the evolution of the toroidal field 
component. 



4 SUMMARY 

In this work, we carried out collapse simulations of magne- 
tised and rotating molecular cloud cores using the Smoothed 
Particle Magnetodynamics (SPMHD) method. We repro- 
duced important features of protostellar outflows, qualita- 
tively in good agreement with results reported in the lit- 
erature. Furthermore, we computed the ratio Btov/B po \ and 
thus give a quantitative proof of the evolution of the toroidal 
part of the magnetic field connected with the outflow. 
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